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Abstract 

Background: Feedback loops, both positive and negative are embedded in the Mitogen Activated Protein Kinase 
(MAPK) cascade. In the three layer MAPK cascade, both feedback loops originate from the terminal layer and their 
sites of action are either of the two upstream layers. Recent studies have shown that the cascade uses coupled 
positive and negative feedback loops in generating oscillations. Two plausible designs of coupled positive and 
negative feedback loops can be elucidated from the literature; in one design the positive feedback precedes the 
negative feedback in the direction of signal flow and vice-versa in another. But it remains unexplored how the two 
designs contribute towards triggering oscillations in MAPK cascade. Thus it is also not known how amplitude, 
frequency, robustness or nature (analogous/digital) of the oscillations would be shaped by these two designs. 

Results: We built two models of MAPK cascade that exhibited oscillations as function of two underlying designs of 
coupled positive and negative feedback loops. Frequency, amplitude and nature (digital/analogous) of oscillations 
were found to be differentially determined by each design. It was observed that the positive feedback emerging 
from an oscillating MAPK cascade and functional in an external signal processing module can trigger oscillations in 
the target module, provided that the target module satisfy certain parametric requirements. The augmentation of 
the two models was done to incorporate the nuclear-cytoplasmic shuttling of cascade components followed by 
induction of a nuclear phosphatase. It revealed that the fate of oscillations in the MAPK cascade is governed by the 
feedback designs. Oscillations were unaffected due to nuclear compartmentalization owing to one design but were 
completely abolished in the other case. 

Conclusion: The MAPK cascade can utilize two distinct designs of coupled positive and negative feedback loops to 
trigger oscillations. The amplitude, frequency and robustness of the oscillations in presence or absence of nuclear 
compartmentalization were differentially determined by two designs of coupled positive and negative feedback 
loops. A positive feedback from an oscillating MAPK cascade was shown to induce oscillations in an external signal 
processing module, uncovering a novel regulatory aspect of MAPK signal processing. 



Background 

Signal transduction pathways such as the Mitogen Acti- 
vated Protein Kinase (MAPK) cascade responds to wide 
range of external stimuli to trigger growth, cell-division 
and proliferation [1,2]. The evolutionarily conserved 
structure of the three layer MAPK cascade consists of 
the MAPKKK (henceforth referred as M3K), MAPKK 
(henceforth referred as M2K) and MAPK (henceforth 
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referred as MK) from yeast to human, which processes 
the incoming signal through a series of covalent modifi- 
cation cycles [1]. M3K is activated upon single phosphor- 
ylation whereas M2K and MK are both activated upon 
double phosphorylation [2-4]. Parallel to the phosphoryl- 
ation by kinases, phosphatases present in the cellular 
volume dephosphorylates the phosphorylated kinases. 
Figure 1 shows the schematics of a three layer MAPK cas- 
cade where each layer of the cascade is dephosphorylated 
by a specific phosphatase. Phosphorylated M3K (M3K*) is 
dephosphorylated by a phosphatase PI, phosphorylated 
forms of M2K (singly phosphorylated M2K* and doubly 
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Figure 1 Schematic representation of the three layer MAPK 
cascade. The cascade is activated by incoming signal represented 
as "Signal". The Signal phosphorylates M3K to M3K* where "*" 
represents phosphorylation. M3K* phosphorylates its downstream 
kinase M2K in two successive steps to M2K** The doubly 
phosphorylated M2K** phosphorylates MK to MK** and MK** is the 
output of the MAPK cascade. The phosphatases PI, P2 and P3 
dephosphorylate the kinases M3K, M2K and MK respectively. The 
rates v1-v10 in a sequential order represents the reaction fluxes. 
Here the solid lines with blunt heads correspond to 
dephosphorylation steps carried out by the phosphatases and 
dashed lines with blunt heads corresponds to phosphorylation of 
the cascade kinases by their respective upstream kinases. 



phosphorylated M2K**) are dephosphorylated by P2 and 
phosphorylated forms of MK (singly phosphorylated MIC 
and doubly phosphorylated MIC*) are dephosphorylated 
by a phosphatase P3. Various feedback loops, both positive 
and negative in nature are abundant in the biological signal 
processing pathways. In a three layer MAPK cascade both 
positive and negative loops are found to be operational [5]. 

Coordinated actions of coupled positive and negative 
feedback loops have been reported earlier for biochemical 
systems with different architectural designs. In cyclin- 
dependent kinase 1 (CDK1) pathway, coupled positive and 
negative feedback loops leads to robust oscillations where 
time periods of oscillations can be changed without com- 
promising the amplitude of oscillations [6]. In another 
study, it was found that during calcium spike regulation, 
positive feedback loops constituting IP3R and RYR and a 
negative feedback loop constituting SERCA ATPases trig- 
gers and regulates the Ca2 ++ oscillations [7]. Similarly the 
cell cycle oscillations are essentially built from coupled 
positive and negative feedback loops between Cdc2 and 
APC system that gives reliable cell cycle oscillations [8]. 



The p53 pathway which is another oscillatory pathway is 
also densely wired in positive and negative feedback loops 
[9]. A recent experimental study shows that a negative 
feedback from the phosphorylated ERK (MK**) to its up- 
stream activator SOS and a coupled positive feedback 
from MK** to M3K* (by phosphorylation mediated dis- 
sociation of M3K* complex from its inhibitor RKIP) 
results in robust system-level oscillations [5], suggesting 
for the first time that the MAPK pathway can employ 
coupled positive and negative feedback loops for generat- 
ing its oscillations. 

In the three-layer MAPK cascade, both positive and 
negative feedback loops emerges from the fully phos- 
phorylated MK (MK**). Feedback loops from MK** act on 
its upstream M2K and M3K layers and alter their phos- 
phorylation according to the nature of the feedback loop. 
A list of feedback loops reported to be operative between 
MK & M2K or MK & M3K are listed in Table 1. It implies 
from the Table 1, that two distinct designs of coupled posi- 
tive and negative feedbacks can potentially exist in the 
three layer MAPK cascade. One design comprises a nega- 
tive feedback from MK** to M3K phosphorylation coupled 
to a positive feedback from MK** to M2K phosphorylation, 
which we named as PN-I design. The other design shows a 
positive feedback from MK** to M3K phosphorylation 



Table 1 Schematics of positive and negative feedback 
loops operational in MAPK cascade 

Feedback interaction of kinases Feedback types References 



Sub module 1 



Positive 



25, 46 




Sub module 2 



Negative 



41,47 




Sub module 3 



Positive 




I MK " H 



Sub module 4 



Negative 



24 
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coupled to negative feedback from MIC* to M2K phos- 
phorylation which is represented as PN-II design. Although 
it is observed that coupling of both positive and negative 
feedback loops can trigger oscillations in the MAPK cas- 
cade [5], potential of both the designs for generating oscilla- 
tions in the MAPK cascade remains to be elucidated. 

Also during long duration signaling, MK and its phos- 
phorylated forms (MK* and MK**), traverses between cyto- 
plasm and nucleus [10-12]. Inside the nucleus, MK** induces 
expression of its phosphatase (MKP-1) that subsequently car- 
ries out MK** dephosphorylation in the nucleus itself [4,10]. 
It is not known how nuclear-cytoplasmic shuttling of the ter- 
minal layer kinase of MAPK cascade and the subsequent 
transcriptional induction of phosphatase such as MKP-1 
would affect the oscillations triggered by PN-I and PN-II. 

Here we built two oscillating models of MAPK cascade 
where oscillations in one model were triggered by PN-I 
and the oscillations in the other model were triggered by 
PN-II. We found that in both the cases, the amplitude, fre- 
quency and nature (digital/analogous) of oscillations were 
uniquely shaped by the coupled positive and negative feed- 
back loops embedded in the cascade. Our simulations 
show that the MAPK cascade embedded in PN-II exhib- 
ited remarkable robustness in generating oscillations with 
identical frequency and amplitude while subjected to a 
wide range of input stimuli, whereas, the cascade embed- 
ded in PN-I was less robust in maintaining its frequency 
and amplitude when subjected to input signal of different 
strengths. We also found that a positive feedback emerging 
from an oscillating MAPK cascade and functional in a dif- 
ferent pathway or signaling module could lead to both sig- 
nal amplification and oscillations in the external module. 
Further we investigated the fate of oscillations of the 
MAPK cascade considering the nuclear and cytoplasmic 
shuttling. Our analysis revealed that the oscillations of the 
MAPK cascade embedded in PN-I were not affected by 
such shuttling of the cascade components and induction 
of its nuclear phosphatase, whereas oscillations triggered 
by PN-II were completely abolished when induction of nu- 
clear phosphatase was considered. Sensitivity analysis for 
small perturbations in parameters of the oscillating models 
were carried out which showed that the organization of 
the feedbacks (PN-I or PN-II) also distinctly determines 
the most sensitive kinetic parameters in the oscillating sys- 
tems. Biological significance of our findings is discussed. 

Methods 

I. Model building 

Information in the signaling cascades such as MAPK 
cascade propagates as a result of phosphorylation- 
dephosphorylation of the kinases in the cascade. Upon ex- 
ternal stimulation, M3K is phosphorylated once and the 
M3K* acts as an enzyme in phosphorylating its down- 
stream kinase M2K. M3K* doubly phosphorylates M2K in 



two single phosphorylation steps [3]. Similarly, the doubly 
phosphorylated M2K (M2K**) phosphorylates MK to MK** 
in two steps. MK** is the output of the MAPK cascade. 
Every layer of the cascade has its individual phosphatase 
[3], which carries out the dephosphorylation process con- 
current to the phosphorylation process. Lets assume an 
external signal 'Signal' triggers the phosphorylation of M3K 
and a cellular phosphatase PI dephosphorylates the phos- 
phorylated M3K* back to its unphosphorylated form. 
The biochemical reaction for phosphorylation process is 
given as 

Signal + M3/Co[Signal.M3/C] Signal + M3IC 

And the biochemical reaction for dephosphorylation 
process is given as 

PI + M3K**> [PIM3IC] — ►PI + M3K 

In the phosphorylation reaction, the 'Signal' could be an 
upstream kinase [13] or other activators that triggers M3K 
phosphorylation [3,10]. The phosphorylation-dephosphor- 
ylation cycles follow in the M2K and MK layers and the 
cascade delivers its final output MK**. 

Under the steady state of production and degradation 
of [Signal.M3K] and [P1.M3K*], flux equations of M3K 
phosphorylation and dephosphorylation (vl and v2 in 
Figure 1) can be given as [14,15] 



vl 



v2 



Kl 

(i+W) 

k 2 .Pl.M3K* 
K2 



Sig.M3iC 

Kl 
i | M3K \ 
1 "•" Kl J 



(i) 

(2) 



Where Sig = k\ .Signal, k x and k 2 are the catalytic rates 
associated with the phosphorylation and dephosphoryla- 
tion processes respectively. Kl and K2 are the Km values 
of the reactions. Phosphorylation-dephosphorylation reac- 
tions for the M2K and MK layer takes place in two steps 
and the equations could be derived accordingly assuming 
steady state conditions [14,15]. 

In MAPK cascade, both positive and negative feedback 
loops emerge from MK** and are functional in either of 
the two upstream layers, M2K and M3K (Table 1). Thus 
the flux equations will be modified in presence of these 
feedback loops. 

Phosphorylation of M3K as shown in equation (1) would 
be modified in presence of a negative feedback loop as, 



Sig.M3/C 
Kl 



(3) 



In equation (3), 'Kl' captures the strength of negative 
feedback of MK** on M3K phosphorylation. The nega- 
tive feedback is assumed as a hyperbolic modifier, which 
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is non competitive in nature and 'nT is the associated 
Hill coefficient [14]. The subscript neg' associated with 
vl in equation (3) represents phosphorylation in pres- 
ence of negative feedback. 

In presence of the positive feedback loop, the flux of 
M3K phosphorylation is modified as 



vl 



(sig-m(i+(^Fr) 



pos 



1 + 



M3K 
Kl 



).(i+erf) 



(4) 



In equation (4), A and Ka are the kinetic constants 
associated with the positive feedback from MIC to the 
M3K layer phosphorylation [16]. The subscript pos' 
associated with vl in equation (4) represents phosphor- 
ylation in presence of positive feedback. In equation (4) 
the exponent nl is the Hill coefficient which indicates 
that the positive feedback is a hyperbolic modifier of the 
M3K phosphorylation [16]. The positive feedback was 
assumed as a hyperbolic modifier in all the model equa- 
tions involving the positive feedback. 

A set of coupled ordinary differential equations cap- 
ture the signal flow in the MAPK cascade which are 
given as 



d[M3IC] 

dt 
d[M2IC] 

dt 

d[M2IC* 

dt 

d[MK*\ 

dt 

d[MK**] 



vl - v2 



v3 - v4 + v5 - v6 



v4 - v5 



vl - v8 + v9 - vlO 



v8 - v9 



The vi, i = 1-10 are the flux equations as given in Table 2 
and also shown schematically in Figure 1. The amount of 
M3K, M2K and MK at any point of time can be calculated 
from the following mass conservation equations. 



[M3K} Tota[ = [M3K] + [M3K*] 
[M2K] Total = [M2K] + [M2K*] + [M2K**] 
[MK] Total = [MK] + [MIC] + [MK**] 

As the total concentration of a kinase is known, M3K, 
M2K and MK can be calculated from the above mass 
conservation equations and the differential equations. 

Models SI and S2 

Based on different types of positive and negative feedback 
loops reported (shown in Table 1), two distinct designs of 



Table 2 Flux of signal flow in cytoplasmic MAPK cascades 'SI' and 'S2' 



Model reactions 



Flux equations in model SI 



Flux equations in model S2 



1] M3K— >M3K* 
2] M3K*— ► M3K 
3] M2K— >M2K* 
4] M2K*->M2K** 
5] M2K**— > M2K* 
6] M2K*— > M2K 
7] MK^ MK* 
8] MK*^ MK** 
9] MK** -> MK* 
10] MK*^ MK 



rpOS " 



Vlneg = 

_ k2., 
v2 = — 

v3 t 

v4, 

v5 = 

v6 = 

v7 -- 

v8 = 

v9 = 

vlO 



(«./W3/(*.^).(l + (^) n: 



(l+^+^)-(l + (^) n3 ) 
(H.M3r.Mfl).(l + (^) n4 ) 

■(i+m n4 ) 



k5.P2.M2K** 



1 K8 T /(7 
k9.P3.M3K* A 



_ (si g .^-) .(i+m n1 
pos ~ (i*w).(i+wsr 



v2 - 



1 1+^ 



neg (i+^r+«-(i+m n3 ) 



v4 nt 



V5 : 



v6~- 



V7 : 



V8 : 



v9-- 



(H./W3K*.^f) 



" (i+^T+^).(i+m" 



''la'' K8 

k8.M2K** MK* 
. K8 
' 1 I MK* I MK 
1 ' K8 'K7 

k9.PZ.MK** 



\/10 : 



The fluxes of signal flow and the values of kinetic parameters used for simulation of SI, S2, S1n and S2n are shown in Table 2. In this table, Ki, i = 1-10 are the Km 
values of the reactions and ki, i = 2-10 are the kcat values of the reactions. The numerical value of T corresponding to Ki and ki represents the reaction number, Kl 
are the kinetic parameters associated with negative feedback. Ka and A are the kinetic constants associated with the positive feedback. The hill coefficient used in 
the equations 1, 3 and 4 are shown as n1, n3 and n4 respectively. 



Sarma and Ghosh BMC Research Notes 2012, 5:287 
http://www.biomedcentral.eom/1756-0500/5/287 



Page 5 of 20 



coupled positive and negative feedback loops emerges. 
Figure 2 shows the two designs of coupled positive and 
negative feedback loops functional in a MAPK cascade, 
namely SI and S2. SI comprises negative feedback from 
MK** to M3K layer coupled to positive feedback from 
MK** to M2K layer. In S2, negative feedback from MK** 
to M2K layer is coupled to positive feedback from MIC* 
to M3K layer. The flux equations of models SI and S2 are 
given in Table 2. All the flux equations corresponding to 
dephosphorylation are identical to each other in both SI 
and S2. Also the flux equations of phosphorylation corre- 
sponding to MK layer are identical in both SI and S2. 
Both SI and S2 were simulated to understand the signifi- 
cance of PN-I and PN-II designs in generating oscillations 
in the MAPK cascade. We studied the characteristic fre- 
quency, amplitude and robustness of the oscillations trig- 
gered by designs, PN-I and PN-II. 

Modification of the models SI and S2 to incorporate 
nuclear-cytoplasmic shuttling 

Nuclear-cytoplasmic shuttling of the MK layer compo- 
nents (MK, MIC and MIC*) of the MAPK cascade takes 
place [10-12] where MIC* triggers various transcription 
factors in the nucleus, aiming to activate target genes 
[10,17]. We updated the models SI and S2 to Sin and S2n 
respectively, to incorporate the nuclear-cytoplasmic trans- 
location of the MK layer components of the cascade. 

In both the modified models, MK** translocate to the 
nucleus and induces its own phosphatase MKP-1 
(named P3-n in the models). The biochemical reactions 
and flux equations corresponding to MK layer s nuclear- 



cytoplasmic shuttling and the transcriptional induction 
of P3-n were adopted from a recent study [10], which is 
given in Table 3. The models Sin and S2n comprise of 
22 flux equations where the first 10 equations in Sin 
and SI are identical to each other which are given in 
Table 2. Similarly the first 10 flux equations of model 
S2n are identical to that of model S2 (Table 2). The add- 
itional equations shown in Table 3 incorporates the nu- 
clear cytoplasmic shuttling of the MK layer components 
MK, MK* and MK**. These also include the equations 
that capture the induction of mRNA of P3-n from the 
target gene triggered by MK** in the nucleus and the 
subsequent biochemical steps that leads to P3-n produc- 
tion. The transcriptionally induced phosphatase P3-n 
dephosphorylates MK* and MK** in the nucleus. The 
differential equations corresponding to the modified sec- 
tion of the model can be found in the Additional file 1: 
model files Sin and S2n. The mass conservation equa- 
tions are identical for SI, S2, Sin and S2n. 

II. Model assumptions 

In substantiation with the previous studies [14,15], it 
was assumed that a steady state in the enzyme-substrate 
complexes is achieved during the signal propagation, for 
all the reactions in both SI and S2. For the sake of sim- 
plicity we assumed that no degradation and production 
of the cascade components (kinases and phosphatases) 
of SI and S2 takes place during the course of signal 
propagation and hence their concentrations remain con- 
stant. However, following experimental guidelines [10], 
the models Sin and S2n were built with certain 



Signal 



MK 
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[M3K* 



|M2K^ 



II 

[ MK* I 
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M2K* 1 
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Figure 2 Schematic representation of the cascades SI and S2. The MAPK cascades embedded in two differential designs of positive and 
negative feedback loops are shown here. (A) An external signal 'Signal' activates the cascade SI by triggering the phosphorylation of M3K. The 
cascade is embedded in the feedback design PN-I. Negative feedback from MK** to M3K layer is shown with a red bar with blunt head and 
positive feedback from MK** to M2K layer is shown with blue arrows. (B) The external signal 'Signal' triggers M3K phosphorylation and the 
cascade S2 embedded in the feedback design PN-II is activated. Positive feedback from MK** to the M3K layer and negative feedback from MK*' 
to the M2K layer are shown as blue arrow and red bars, respectively. 



Sarma and Ghosh BMC Research Notes 2012, 5:287 
http://www.biomedcentral.eom/1756-0500/5/287 



Page 6 of 20 



Table 3 After addition of the nuclear components, 'SI' and 'S2' were renamed as 'Sin' and 'S2n' respectively 



Reaction number in models Sin and S2n 


Reaction 


Flux equation 


Parameter value 


11] 


MK** <-> MK**-n 


k11 f .MK**- k11 b .MK**-n 


k11 f = 10.34 sec" 1 
k1 1 b = 2.86 sec" 1 


12] 


P3-n gene PreP3mRNA 


l/12.(/W/C*-n) n12 
(/(12) n12 +(/M/(**-n) n12 


VI 2 = 29.24 nmol/sec 
K12 = 169 nmol/ml 
n12 = 3.97 


13] 


PreP3mRNA^P3mRNA 


k13. PreP3mRNA 


kl 3 = 0.022 sec" 1 


14] 


P3mRNA -> (|) 


k14.P3mRNA 


kl 4 = 0.0078 sec" 1 


15] 


P3mRNA -> P3-c 


k15.P3mRNA 


kl 5 = 0.001 2 sec" 1 


16] 


P3-c <-> P3-n 


k16f.P3-c- k16 b .P3-n 


k16f= 22.56 sec" 1 
k16 b = 15.4 sec" 1 


17] 


P3-c -> (|) 


k17.P3-c 


k1 7 = 0.00025 sec" 1 


18] 


P3-n -> (|) 


k18.P3-n 


kl 8 = 0.00025 sec" 1 


19] 


MK^ MK-n 


k19 f .MK- k19 b .MK-n 


k19 f = 10.34 sec" 1 
k19 b = 2.86 sec" 1 


ZUJ 


IVIix > I VI IX II 


KZUf.lVIIX KZUf.lVIIX II 


ls~)O r — 1 n cor" 1 
KZUf — 1 U. ScL 


21] 


MK**-n -> MK*-n 


*21.(P3-n).^> 

11 K21 ' /C22 


k20 b = 2.86 sec" 1 

k21 =0.68 sec" 1 

K21 =10300 nmol/ml 


22] 


MK*-n -> MK-n 


*22.(P3-n).^ 

-i , MK* — n , /MK**-n 
11 K22 1 K21 


K22 = 87 nmol/ml 
k22 = 0.31 sec" 1 
K21 =10300 nmol/ml 
K22 = 87 nmol/ml 



Table shows the reactions and flux equations in both S1n and S2n are shown in Table 3. Reactions 1-10 in Sin and S2n were identical to that given in Table 2 
for S1 and S2 respectively. Reactions and flux equations shown in this table are common for both S1n and S2n. Reactions 21 and 22 were built from the original 
mass action reactions ([10], by assuming steady state condition of the enzyme-substrate complexes, which were used to calculate Km value for reactions 21 and 
22. In the reactions, ki f and ki b where i = 1 1, 16, 19, 20 are the forward and backward reaction rates of the respective reactions. The parameter values ki, where 
i = 13, 14, 15, 17, 18, 21, 22 are the catalytic rates of the respective reactions. Ki, i = 12, 21, 22 are the Km values of the respective reactions. V12 is the maximum 
velocity of the reaction 1 2. 



degradation and phosphatase (P3-n) production steps, as 
shown in Table 3. In models SI and S2 we also assumed 
that each layer of the cascade is phosphorylated by one 
phosphatase specific to each layer [3,14], except, in the 
models Sin and S2n, where dephosphorylation of the 
third layer MK was carried out by two phosphatases, P3 
and transcriptionally induced P3-n. The model presented 
here represents a three layer MAPK cascade that is evolu- 
tionarily conserved from yeast to mammal [1]. Although 
differences in the rewiring of the kinases-phosphatases 
interaction are observed in some eukaryotic systems 
[13-15,18], the kinases-phosphatases interaction shown 
here represents the most generalized structure of the 
cascade known till now [1-3]. The simplifications also 
included ignoring various intra modular crosstalks which 
involve MAPK cascade and other signaling modules [19]. 
While building the flux equations for positive and negative 
feedback loops we assumed that both the feedback types 
are hyperbolic modifiers, which is in corroboration with 
earlier studies [14,16]. 



III. Model parameters and concentrations 

The kcat and Km values for SI, S2, Sin and S2n were 
chosen in biochemically observed ranges [3,14-16,18]. 
Additional file 2: Table SI describes the reactions cap- 
turing signal flow in the three layer MAPK cascade and 
their kinetic parameter values, which are common in all 
the four models SI, S2, Sin and S2n. Additional file 2: 
Table S2 describes the concentration of kinases and 
phosphatases used in SI, S2, Sin and S2n. Table 3 shows 
the additional reaction parameters corresponding to the 
modified fraction of the models Sin and S2n. Para- 
meters for the additional reactions in the model Sin and 
S2n were adopted from a recent study [10]. 

IV. Sensitivity analysis for small perturbations in the 
model parameters 

Sensitivity studies reveal the relative importance of kin- 
etic parameters associated with the model. We per- 
formed sensitivity analysis of all the four models by 
applying small perturbations to the kinetic parameters of 
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the models and measuring the sensitivity of MIC* in 
each of the model to such perturbations. 

Mathematically, the sensitivity coefficients are the first 
order derivatives of model outputs with respect to the 
model parameters; Sg = , where Oi is the i th model 

output and pj is the j th model parameter [20,21]. Sy is 
the sensitivity coefficient which yields sensitivity of Oi 
with respect to the perturbation in parameter p^. We 
have calculated the sensitivity coefficient Sy using the 
software SBML-SAT that implements the centered dif- 
ference assumption for calculating [21]. When a par- 
ameter pj is subjected to a small perturbation (Apj) in its 
reference value (reference value is the unperturbed 
value), the sensitivity coefficient is calculated as 

q _ Oi(pj+Apj)-O i (p j -Apj) 
" ~ 2A PJ 

Upon normalization, the sensitivity coefficient is 

Oi (pj+Apj ) - O t (pj -A P j ) 

given as: Sy normalized = 2 %. - 

In the above equation, we calculated Sy with Apj = O.OOPpj 
for any perturbed parameter py The variation of Apj in the 
range of 0.0001*pj-0.1*pj didn't alter Sy. The perturbations 
were applied locally, which means parameters were 
perturbed one at a time and for each of the parameter s 
perturbation on the output MIC* of the models was 
calculated. 

V. Software used and model simulations 

For performing the simulations SBML models were ini- 
tially constructed using Complex pathway simulator 
(COpasi) [22]. The time course simulations were carried 
out in COpasi. Sensitivity analysis was performed using 
SBML_SAT, a MATLAB toolbox for sensitivity analysis 
[21]. Bifurcation analysis to inspect oscillation in S2n 
was carried out using Bifurcation Discovery tool [23]. 
The model files are given as additional files. 

Results 

We constructed two models SI and S2 of the MAPK 
cascade, one embedded in PN-I and the other embedded 
in PN-II respectively, such that oscillations in both the 
models were triggered by coupled positive and negative 
feedback loops. We investigated the fate of MAPK oscil- 
lations in SI and S2, when signal strength was varied in 
wide ranges. Our simulations also revealed that MAPK 
cascade can utilize its positive feedback to trigger oscilla- 
tions in an external signal processing module. Next we 
examined the fate of oscillations triggered by PN-I and 
PN-II when nuclear-cytoplasmic shuttling of the com- 
ponents of terminal layer MK of the MAPK cascade 
takes place followed by the induction of a nuclear 



phosphatase by MK**. Results show that oscillations 
triggered by PN-II exists only in the cytoplasm and in- 
duction of the P3-n completely abolished the oscilla- 
tions, whereas oscillations triggered by PN-I are not 
affected by the nuclear translocation of MK layer and 
subsequent induction of nuclear phosphatase. Various 
in-silico knock out studies were carried out to elucidate 
the importance of cytoplasmic and nuclear phosphatases 
in both SI and S2. Also, when the parameters of SI, 
Sin, S2 and S2n were subjected to small perturbations, 
we found that PN-I and PN-II differentially regulates the 
cascades' output sensitivity to these perturbations. 

Oscillations in models SI and S2 

Previous studies show that negative feedback from MK** 
to M3K layer (Sub module 2, Table 1) [14], or negative 
feedbacks from MK** to M2K layer (Sub module 4, 
Table 1) [24], triggers sustained oscillations in the 
MAPK cascade. Positive feedback from MK** to M3K 
phosphorylation (Sub module 1 Table 1) results in all- 
or-none behavior in production of MK** [4,25,26]. Posi- 
tive feedback from MK** to M2K phosphorylation step 
(Sub module 3, Table 1) was found to facilitate propaga- 
tion of long range phosphorylation waves of MK** in the 
developing neurons [16]. Earlier computational investi- 
gations revealed that a negative feedback from MK** to 
M3K layer is a prerequisite in triggering MAPK oscilla- 
tions [14], but later it was found that for certain param- 
eter combinations, the three layer MAPK cascade can 
trigger its oscillations in absence of the explicit negative 
feedback loop from MK** to M3K [27]. 

But a recent experiment exposed that MAPK oscilla- 
tions are triggered by coupled positive and negative feed- 
back loops [5]. This experimental finding necessitated an 
investigation on the significance of differential designs of 
coupled positive and negative feedback loops that can 
plausibly trigger oscillations in the cascade and the char- 
acteristics of oscillations triggered by each of the design. 
The MAPK cascades embedded in the two designs of 
coupled positive and negative feedback loops, PN-I and 
PN-II are shown in Figure 2A and 2B. 

Upon simulation of models SI and S2 without any 
feedback loops, maximum amplitude phosphorylation of 
the output (MK**) was attained (Figure 3A and 3D). 
When both the models were simulated in presence of 
only negative feedback loops, MK** amplitude was 
inhibited (Figure 3B and Figure 3E). This occurs as a re- 
sult of negative feedback mediated suppression of M3K 
layer phosphorylation in SI and M2K layer phosphoryl- 
ation in S2 respectively. We show that the employed 
negative feedbacks lead to inhibition of MK** amplitude 
in both SI and S2, to demonstrate that the models con- 
sidered in our studies do not oscillate only in the pres- 
ence of the negative feedback loop like the earlier 
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Figure 3 Dynamics of MAPK cascade embedded in PN-I and PN-II. Both models SI and S2 were simulated without any feedbacks, with 
negative feedback alone and with coupled positive and negative feedbacks. Without any feedback loop, both SI and S2 shows full amplitude 
phosphorylation of MK**, M2K** and M3K* (A and D). The negative feedback from MK** to M3K phosphorylation in SI [B] and from MK** to M2K 
phosphorylation in S2 (E) resulted in inhibition of MK** amplitude. The coupled feedback design PN-I in SI (C) and PN-II in S2 triggered sustained 
oscillations (C & F). 



reports [14,24]. Next when positive feedbacks were 
introduced in the models, both SI and S2 exhibited sus- 
tained oscillations (Figure 3C and 3F), demonstrating 
that the MAPK cascades considered for our study oscil- 
late only in the presence of coupled positive and nega- 
tive feedback loops. 

Oscillations in SI 

Introduction of the positive feedback loop from MK** to 
M2K layer in the cascade with negative feedback from 
MIC* to M3K layer first resulted in enhancement of the 
amplitude of M2K** followed by enhancement in MK** 
amplitude. Since both positive and negative feedbacks 
emerges from MK**, enhanced MK** amplitude results 
in stronger inhibition in the M3K* layer and stronger ac- 
tivation in the M2K layer. However as M3K lies up- 
stream to M2K, decrease in M3K* concentration beyond 
a certain threshold results in attenuation of M2K layer 
phosphorylation, even in the presence of the positive 
feedback loop. With inhibition of M2K** amplitude, 
phosphorylation of MK layer gets inhibited. With 



decrease in MK layer phosphorylation, attenuation of 
the strengths of both positive and negative feedback 
loops follow. As MK** amplitude reaches its lowest amp- 
litude, one cycle of oscillation is completed (Figure 3C). 
As the input signal is available for M3K phosphorylation, 
M3K* starts building up in absence of the negative feed- 
back and the next cycle of oscillation is triggered. The 
process continues until the external signal is available to 
phosphorylate M3K. Coupling of inhibitory and activat- 
ing effects of the PN-I, triggered oscillations (Figure 3C) 
in all the three kinases of the MAPK cascade SI. 

Oscillations in S2 

Oscillations in S2 emerged due to positive feedback 
mediated enhancement of M3K* amplitude coupled to 
the negative feedback mediated inhibition of M2K**. 
Upon stimulation of the cascade by external signal, posi- 
tive feedback from MK** to M3K enhanced the M3K* 
amplitude. This subsequently enhances M2K layer phos- 
phorylation (in presence of the negative feedback from 
MK** to the M2K layer), ultimately resulting in 
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amplification of MIC* amplitude. Amplified MIC* subse- 
quently enhances the strengths of both positive and 
negative feedback loops. When MIC* reaches its max- 
imum phosphorylation amplitude (Figure 3F), negative 
feedback mediated inhibition of M2K layer phosphoryl- 
ation surmounts the positive feedback mediated en- 
hancement of M2K layer phosphorylation by M3IC. 
With progressive attenuation of M2IC* amplitude, MK 
layer phosphorylation gets inhibited until it reaches its 
lowest phosphorylation amplitude (Figure 3F). The 
whole process completes one cycle of oscillation. The 
next cycle of oscillation starts when the external signal 
triggers phosphorylation of M3K in absence of the nega- 
tive feedback from MK**. It could be noted that the 
negative feedback in S2 inhibits MIC* production in two 
ways, firstly by directly inhibiting the M2IC* amplitude 
and secondly by indirectly inhibiting the M2IC* by at- 
tenuating the strength of positive feedback loop from 
MIC* to the M3K layer. The study additionally uncov- 
ered that positive feedback not only enhanced M3IC 
amplitude but it also triggered oscillations in M3IC 
(Figure 3F). 

Nature of oscillations in SI and S2 

In SI, where the incoming signal encounters the nega- 
tive feedback first and then the positive feedback, output 
oscillations (MK**) are digital in nature (Figure 3C). In 
S2, the signal encounters positive feedback first followed 
by its encounter with the negative feedback, which 
resulted in sinusoidal oscillations (Figure 3F). In the 
MAPK cascade, it is known that positive feedback stabi- 
lizes [25,28] and negative feedback destabilizes [14,24] 
the output (MK**) amplitude. Here we showed that the 
interplay between such stabilizing and destabilizing ef- 
fect differentially determines the nature of oscillations 
which ultimately depends on the designs of coupled 
feedback loops. The digital oscillations in SI exhibited 
sharp switch like characteristics of a positive feedback 
[29] in the rise and fall of the phosphorylation waves 
(Figure 3C) and the analogous oscillations in S2 exhib- 
ited characteristics of a negative feedback mediated 
oscillations observed earlier [11]. The study suggests that 
output characteristics of an oscillating MAPK cascade is 
based on the feedback type encountered by the incoming 
signal at the M2K layer. 

Next we examined how oscillations in the MAPK cas- 
cade embedded in PN-I and PN-II are affected when 
both SI and S2 are activated by input signal of different 
strengths. 

Oscillations in SI and S2 subjected to a wide range of 
input stimuli 

Signal strength varies widely in the in-vivo conditions. 
The strength of the incoming signal is governed by the 



concentration of the signal as well as the proximity of 
the signal source to the target receptor that activates a 
signaling pathway [3,5,30,31]. However biological sys- 
tems are built to maintain their output characteristics in 
the face of perturbations [32]. Thus we examined the 
relative robustness of SI and S2 in triggering their char- 
acteristic oscillations when both the systems were sub- 
jected to a spectrum of input signals. 

I. Model SI 

Figure 4A shows the oscillation characteristics of SI sub- 
jected to a range of input signals (varied between 2.5 nM- 
50 nM; 2.5 nM is the minimum signal strength required 
for triggering the oscillations). At a low signal strength, 
(shown as 'Sig' in Figure 4A) MK** oscillations with 
maximum amplitude were achieved. With increase in 
signal strength, the effect of negative feedback mediated 
suppression of M3K phosphorylation was diluted and 
beyond a certain strength of the input signal (Sig > 50 
nM), the negative feedback can no longer suppress 
M2K layer phosphorylation by inhibiting M3K phos- 
phorylation. Thus beyond a certain strength of input 
signal (Sig > 50 nM), coupled effect of the strong input 
signal and the positive feedback from MK** to M2K 
layer resulted in a steady non-oscillatory phosphoryl- 
ation of M2IC* and MK** (data not shown). However if 
the signal was applied in the range provided above, sus- 
tained oscillations could be achieved in the cascades 
output phosphorylation (Figure 4A). With increase in 
signal strength (in the range 2.5-50 nM), oscillation 
amplitudes were conserved, but the frequency of oscil- 
lations decreased with increasing strengths. Thus a 
MAPK cascade embedded in PN-I can exhibit con- 
served amplitude oscillations whose frequencies would 
be decided by the strengths of the incoming signal. 

II. Model S2 

The model S2 was subjected to signals of variable 
strengths. Beyond a certain threshold (Sig ~ 5 nM) that 
triggered oscillations in the cascade, oscillations were 
observed for signals of any given strength (we tested the 
Sig range in 5 nM - 5000000 nM) of incoming signal. 
Figure 4B shows MIC* oscillations in S2 for the signal 
strength 5-500 nM. S2 also exhibited sustained oscilla- 
tions with equal frequency and amplitude for all the 
strengths of applied signal above the threshold strength. 
The causality behind emergence of such robust oscilla- 
tions could emerge from the design of the coupled feed- 
back loops. In S2, positive feedback enhances M3IC 
amplitude and thus for a relatively smaller signal dose 
M3IC reaches its maximum amplitude and saturates. 
Hence when the signal strength is increased further, no 
additional changes will be observed in the M3K layer. 
Since the strengths of the feedback loops becomes 



Sarma and Ghosh BMC Research Notes 2012, 5:287 
http://www.biomedcentral.eom/1756-0500/5/287 



Page 10 of 20 




Time (seconds) 



Figure 4 Sustained oscillations in SI and S2 as function of signal dose. Strength of the signal activating the MAPK cascade 'Sig' was varied 
and the range of 'Sig' in which both SI and S2 would exhibit oscillations was studied. (A) After a certain threshold value of 'Sig', increase in 'Sig' 
resulted in decrease in the frequency of MK** oscillations in SI, whereas the maximum amplitude of oscillations remained unaltered. When the 
range of 'Sig' triggering such oscillations was crossed beyond 50nM, the oscillations ceased to exist (data not shown). (B) In S2, beyond a 
threshold value of Sig (~ 5 Nm) required to trigger oscillations, almost any value of 'Sig' triggered oscillations with conserved frequency and 
amplitude. 



unresponsive to the further increases in signal strength, 
MK** oscillations with robustly conserved amplitude 
and frequency could be generated for a very wide range 
of input signals. 

As shown earlier for the system S2, positive feedback 
led to oscillations in the M3K* amplitude in addition to 
the amplification in its phosphorylation (Figure 3F). We 
next investigated whether the positive feedback compo- 
nent of S2 (and also SI) is capable of transferring oscilla- 
tions to external signal transduction modules in general. 

Positive feedback transfers oscillations from an oscillating 
MAPK cascade to other signaling modules 

Results shown in Figure 3F opens up a possibility that 
positive feedback loop emerging from an oscillating 
MAPK cascade could trigger oscillations in its place of 
action in addition to the signal amplification in the tar- 
get site. Experimentally such positive feedback loop is 
observed from the output MK** (from p38MAPK cas- 
cade) to the p53 phosphorylation step [9]. Similarly posi- 
tive feedback from the output MK** (from ERK cascade) 
leads to modification of Lck kinase as observed in the 



T lymphocytes [33]. We investigated how the positive 
feedback from oscillating MAPK cascades such as SI or 
S2 would affect the phosphorylation in an external signal 
transduction module, by building a hypothetical 
phosphorylation-dephosphorylation cycle with a kinase 
X and its phosphorylated form X-P. The model used for 
simulation of the positive feedback from S2 to X is pro- 
vided as an additional SBML model file. 

We built a model where MK** of system S2 provides a 
positive feedback to the phosphorylation of a kinase X 
(X is a hypothetical kinase phosphorylated by a different 
signal). Kinase X was assumed to be activated by phos- 
phorylation like most of the kinases in the signaling net- 
works. Also we assumed that a cellular phosphatase 
dephosphorylates phosphorylated X (X-P) back to its 
unphosphorylated form. This simple one step covalent 
modification cycle represents the most fundamental 
module of signal transduction [34] and is a building 
block of nearly all the signal processing modules [35,36]. 

We introduced the positive feedback loop from S2 to 
phosphorylation step of X. Simulations show that the 
positive feedback transfers oscillations from S2 to X 
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where the extent of oscillations in X-P (or X) was gov- 
erned by the relative rates of phosphorylation and 
dephopshorylation in the X module. Figure 5A shows 
the dynamics of X-P phosphorylation in presence and in 
absence of the positive feedback loops, when phosphor- 
ylation rate of X (X phos ) is equal to dephopshorylation 
rate of X-P (X dephos ). The positive feedback transferred 
oscillatory information from S2 to X-P together with 
triggering amplification in X-P amplitude. When the 
X phos is higher than the X dep hos (Figure 5B) oscillatory 
phosphorylation of X was diluted but the amplification 
of X-P caused by the positive feedback remained un- 
affected. On the contrary when X phos < X dephos 
(Figure 5C), X-P exhibited oscillations with much wider 
differences in the maximum and minimum amplitudes 
of its oscillations. For significantly lower values of X phos 
(Xphos < < X dephos ), phosphorylation of X oscillated be- 
tween its lowest (zero) to its maximum phosphorylation 
amplitude (Figure 5D). We also investigated the effect of 
positive feedback emerging from MAPK cascade SI and 
functional in the phosphorylation step of the module X. 
Here too oscillations from SI to the module X were 
transferred as function of relative values of X phos and 



X d ephos with maximum amplitude oscillations in the X 
module triggered when X phos < < X dep hos (data not 
shown) This study exposes a novel cellular strategy 
where cells can control the effects of a positive feedback 
loop emerging from a MAPK cascade such as SI or S2 
and operational on different target sites. We revealed 
how adjustment of phosphorylation and dephosphoryla- 
tion rates in the target modules would regulate the ex- 
tent of oscillations in them. 

Next we investigated the fate of oscillations triggered 
by PN-I and PN-II when nuclear cytoplasmic shuttling 
of the MK layer takes place. The analysis was performed 
to investigate the fate of oscillations triggered by PN-I 
and PN-II when the oscillations in the cascade output 
(MK**) are triggered in the cytoplasm but its nuclear 
translocation takes place subsequently. 

Fate of MAPK oscillations in SI and S2 upon nuclear 
translocation of the MK layer followed by induction of its 
own nuclear phosphatase 

It was observed experimentally that upon prolonged sig- 
naling, nuclear cytoplasmic shuttling of the MK layer of 
the MAPK cascade takes place [10,11]. Activation of the 
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Figure 5 Oscillations in an external signaling module of kinase X by the positive feedback of the MAPK cascade S2. (A) Dynamics of X-P 
with and without the positive feedback loop when rate of phosphorylation of X (X phos ) = rate of dephopshorylation of X (X dephos ). Here 
Xphos = Xdephhos = 0.1 • (B) Dynamics of phosphorylated X (X-P) when X phos > X dephos The X phos = 0.5 and X dephos = 0.1 . (C) Dynamics of X-P when 
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MAPK cascade is followed by nuclear translocation of 
its output MK** where it induces various transcription 
factors including its own phosphatase. It is known that 
upon nuclear translocation, the doubly phosphorylated 
ERIC* (MIC* in our models) induces its nuclear phos- 
phatase MKP-1 [10,37]. The phosphatase MKP-1 is nu- 
clear specific; thus it dephosphorylates MIC* only in the 
nucleus. Hence for the long duration MAPK signaling, 
where induction of the MAPK phosphatase MKP-1 takes 
place [4,37], the phosphorylated MIC* is dephosphory- 
lated in the cytoplasm by P3 and also in the nucleus by 
MKP-1 (P3-n in Figure 6). Here, first we investigated the 
sustainability of oscillations upon nuclear cytoplasmic 
shuttling of the MK layer components and subsequently 
studied the roles of P3 and induced P3-n in determining 



the oscillatory fate of MIC* and its nuclear component 
MIC*-n ( Figure 6). For that purpose the existing models 
SI and S2 were modified to incorporate nuclear trans- 
location of the MK layer and induction of P3-n by 
MIC*-n (Figure 6). The modified models had 22 bio- 
chemical reactions each (SI modified to Sin and S2 
modified to S2n respectively), with the first 10 reactions 
in Sin and S2n being identical to SI and S2 respectively, 
which are shown in Table 2. The 11 additional reactions 
in Sin and S2n (Shown in Table 3) captured shuttling of 
MK, MK* and MK** between cytoplasm and nucleus, 
P3-n induction steps and dephosphorylation of MK**-n 
and MIC-n in the nucleus by P3-n. Mechanistic and 
parametric details for nuclear-cytoplasmic shuttling of 
MK layer components and transcriptional induction of 
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Figure 6 Schematics of compartmentalized MAPK signaling and transcriptional induction of the phosphatase P3-n. The MK, MK* and 

MK** comprises of the cytoplasmic fraction which translocate to the nucleus and their nuclear counterparts are MK-n, MK*-n and MK**-n. The 
cytoplasmic and nuclear components together constitute the total MK. MK**-n transcriptionally induces mRNA of P3-n which is shown as 
PreP3mRNA that finally leads to translation of P3-n in the cytoplasm (shown as P3-c). Upon nuclear translocation P3-c becomes P3-n and 
dephosphorylates MK**-n and MK*-n. Phosphorylation, dephosphorylation, transcription, translation, and nuclear-cytoplasmic shuffling are 
represented with distinct arrowheads/types and lines with blunt heads. 
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P3-n were taken from a recent study on the mammalian 
MAPK cascade [10]. 

I. Oscillations in Sin 

Nuclear compartmentalization of the MK layer and tran- 
scriptional induction of P3-n didn't affect the oscillations 
in Sin and it exhibited MK** oscillations with near iden- 
tical frequencies as observed in SI. Nonetheless, the 
amplitude of cytoplasmic MK** decreased and major 
fraction of phosphorylated MK** resided in the nucleus 
(MK**-n in Figure 7A). Next we checked the roles of P3 
and P3-n in deciding the oscillatory fate of MK** and 
MK**-n. P3 concentration was made 0 (biologically cor- 
responds to a P3 knockout system) and the system was 
simulated. Figure 7B shows the results for P3 = 0, when 
dephosphorylation of MK**-n and MK*-n was carried 
out by P3-n. The simulations show that the frequency 
and amplitude of MK** and MK**-n were not altered 
when P3 is absent in the system and dephosphorylation 
of MK layer is carried out only in the nucleus (Compare 
Figure 7 A and 7B). In the subsequent analysis, we 
stopped transcriptional induction of P3-n after 600 sec- 
onds (randomly chosen) of simulation, with P3 = 0 as an 
initial condition before simulation (In the model, the 
transcriptional induction was stopped by equating the 



rate of P3mRNA production to zero). We observed 
MK** and MK**-n oscillations for significantly long time 
after the transcription was stopped (Figure 7C), and only 
after P3-n concentration goes down a certain limit (due 
to its degradation), oscillations in MK** and MK**-n 
were abolished (Figure 7C). However oscillations could 
be triggered back to the system (Figure 7D) when P3 
concentration was reverted from 0 back to its reference 
value (500 nM, Additional file 3: Table S2) after P3-n 
concentration goes below a value (due to its degrad- 
ation) that is required to maintain sustained oscillations. 
The simulations thus show that MAPK cascade with 
architectural design such as Sin can exhibit oscillations 
in presence of either of the nuclear or cytoplasmic phos- 
phatase. It could be noted that presence of both phos- 
phatases didn't impart any change in the frequencies 
and amplitudes of MK** and MK**-n. 

II. Oscillations in S2n 

Simulations were carried out in S2n after incorporation 
of the transcriptional components (nuclear cytoplasmic 
shuttling and induction of P3-n) in the MAPK cascade. 
Similar to the model Sin, the model S2n was also built 
upon the existing model S2. Similar to Sin, the 
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Figure 7 Oscillations in Sin. Oscillations in the MAPK cascade with underlying feedback design PN-I are not altered or abolished by its nuclear 
translocation and induction of the nuclear phosphatase P3-n (A) Oscillations in MK** and MK**-n when P3 dephosphorylated MK** in the 
cytoplasm and P3-n dephosphorylated MK**-n in the nucleus. (B) Oscillations of MK** and MK**-n when P3 was knocked out from the system, 
implying P3 = 0 throughout the simulations but when P3-n was allowed to be produced transcriptionally. (C) Oscillations in MK** and MK**-n, 
when P3 = 0 throughout the simulations, and P3-n production was stopped at time = 600 seconds. (D) P3 is knocked out initially; P3-n 
production was stopped at 600 seconds, which was followed by reverting the P3 concentration back to its reference value (P3 = 500 nM), at 
10000 seconds. 
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parameters for transcriptional processes were kept iden- 
tical to the experimentally reported values [10]. 

Dynamics of MIC* and MK**-n phosphorylation are 
shown in Figure 8A. The simulations show that when 
MK**-n was used to induce its own phosphatase P3-n, no 
oscillations where observed in the system. When P3 = 0, 
amplitudes of MIC* and MK**-n (Figure 8B), didn't differ 
from the condition when P3 = 500 Nm (Figure 8A). How- 
ever when only the nuclear-cytoplasmic shuttling of MK 
layer components was considered (without P3-n induction 
and keeping P3 = 500 nM), the system exhibited its char- 
acteristic oscillations (Additional file 3: Figure SI). This 
implies that oscillations in S2n were not abolished due to 
nuclear-cytoplasmic shuttling of the MK layer compo- 
nents, but due to the transcriptional induction of P3-n. 
For P3 =0 as an initial condition, followed by inhibition of 
P3-n at 600 seconds, the oscillations in MIC* were not 
observed for any value of P3-n concentration (Figure 8C). 
Next, in a P3 knocked out system, P3-n was made zero at 
time = 600 seconds, followed by reverting P3 concentra- 
tion back to 500 nM after time = 10000 seconds (after 
P3-n was significantly degraded). We found that after 
P3-n concentration becomes significantly low (after time 
= 10000 seconds) reverting P3 back to its reference value 
triggered sustained oscillations in both MK** and MK**-n 
(Figure 8D). Introduction of P3 in presence of higher con- 
centrations of P3-n (time < = 10000) didn't trigger oscilla- 
tions in MK** and MK**-n (Additional file 3: Figure S2). 



We also searched the parameter space (in the range of 
0.01-100 times their reference values as given in Table 3) 
of model S2n for combinations of parameters that could 
possibly trigger sustained oscillations in S2n. The para- 
meters were varied using Bifurcation discovery tool [23] 
where we searched specific combinations of parameters 
that could trigger oscillations in S2n in presence of both 
P3 and P3-n. The analysis provided a parameter set that 
triggered transient oscillations (data not shown), but to 
trigger such oscillations, values of several of the para- 
meters were largely shifted from their experimentally 
observed values. Thus applying changes in those param- 
eter values would perhaps not represent the realistic sce- 
nario anymore and we restricted ourselves from 
applying such changes in S2n. Our analysis thus suggests 
that in a MAPK cascade embedded in feedback design 
such as PN-II, sustained oscillations could only be trig- 
gered in absence of its nuclear phosphatase P3-n. 

PN-I and PN-II differentially shapes the MAPK cascades' 
output sensitivity to small perturbations in parameter 
values 

In signaling networks with multiple parameters, perturb- 
ation in only a few parameters pivotally decides the out- 
put fate of the systems and changes in majority of the 
parameters doesn't alter the output characteristics [38]. 
Knowledge of the crucial and less-crucial parameter 
values improves the understanding on the regulatory 



800 




5600 




•2400 




£200 

0 
c 


V 


0 






2000 



6,000 10,000 15,000 
Time (seconds) 



20,000 




6000 10,000 15,000 
Time (seconds) 



20,000 



800 



E600 



□ 400 



=?200 




0 2000 



6000 10,000 15,000 
Time (seconds) 



20,000 



0 2000 



6000 10000 15000 
Time (seconds) 



20,000 



Figure 8 Oscillations in S2n. Oscillations in the MAPK cascade with underlying feedback design PN-II were terminated due to induction of its 
nuclear phosphatase P3-n. (A) Oscillations in MK** and MK**-n were abolished when P3 dephosphorylated MK** in the cytoplasm and P3-n 
dephosphorylated MK**-n in the nucleus. (B) Oscillations were not observed when P3 was knocked out (P3 = 0) and when only P3-n 
dephosphorylated MK**-n and MK*-n in the nucleus (C) For P3 = 0 throughout the simulations, and P3-n production stopped at time = 600 
seconds, oscillations were not observed for any value of P3-n concentration. (D) When P3 was set to zero as initial condition and P3-n 
production was stopped at 600 seconds, followed by reverting P3 concentration back to its reference value (P3 = 500 nM) at 10000 seconds, 
sustained oscillations in both MK** and MK**-n were observed. 
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principles and helps in finding suitable drug targets 
[39,40]. We subjected the kinetic parameters of SI, S2, 
Sin and S2n to small perturbations and the sensitivities 
of the outputs MIC* (in SI and S2) and MK**-n (in Sin 
and S2n) were calculated. Thus a model parameter 'p' 
was subjected to perturbation Ap where Ap= O.OOPp. 
Such small perturbations in the parameter values didn't 
affect the sustained nature of oscillations, but revealed 
the relative sensitivity of the output to the perturbations. 

Figure 9A and 9B shows the sensitivity of MK** to small 
perturbations in their model parameters. MK** in the 
MAPK cascade embedded in PN-I and PN-II was found 
to exhibit different sensitivity profiles. In the Figure 9A 
and 9B, only the most sensitive parameters are shown 
with their respective names. In SI (Figure 9A), MIC* is 
most sensitive to the perturbations in the strength of the 
incoming signal (Sig) and the dephopshorylation rate (k2) 
of M3K**. In S2 (Figure 9B), MIC* is most sensitive to 
perturbations in rates of dephosphorylation in the MK 
layer (k9 and klO in Figure 9B). The models Sin and S2n 
were also subjected to small perturbations like in SI and 
S2 (Figure 9C and 9D). The sensitivity profile of MK**-n 
in Sin was similar to MIC* in SI with MK**-n being most 
sensitive to changes in signal strength and the 



dephosphorylation rate of M3IC (Figure 9C). MK**-n in 
S2n exhibited relatively higher sensitivities to the para- 
meters involved in the shuttling of MK layer components 
specifically the shuttling rate of MK**-n (Figure 9D). 

The differential sensitivity profile of MIC* in the two 
models could be mechanistically understood as follows. 
The MAPK cascade being a ultrasensitive cascade and 
signal amplifier [3,25], any small changes in the input 
layer gets amplified as it propagates downstream and 
results in significantly larger changes in the output of 
the system. Usually negative feedback is a noise suppres- 
sor and small fluctuations in the values of signal/para- 
meters are filtered by the negative feedback [41,42]. But 
as the positive feedbacks are coupled to the system as 
well they further amplify the effect of small changes/per- 
turbations, and subsequently alter the phosphorylation 
of the MIC*. Thus in SI and Sin (Figure 9A and 9C), 
changes in the M3K layer due to small fluctuations in 
the parameter values were amplified at the M2K layer 
owing to the positive feedback. Thus coupling of the ef- 
fect of the positive feedback together with the MAPK 
cascades inherent ability for signal amplification (due to 
multisite phosphorylation and multi layer organization 
[3,4]) resulted in maximum sensitivity of MIC* to small 
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Figure 9 Sensitivity of MK** in SI and S2 and MK**-n in Sin and S2n to small perturbations in the system parameters. The kinetic 
parameters and signal strength 'Sig' of both SI and S2 were subjected to small perturbations and sensitivity of MK** (in SI and S2) and MK**-n 
(in Sin and S2n) to such perturbations were calculated. (A) Response of MK** of the system S1 to the perturbations in kinetic parameters. 
(B) Response of MK** of the system S2 to the perturbations in kinetic parameters. (C) Response of MK**-n of the system Sin to the perturbations 
in kinetic parameters and (D) Response of MK**-n of the system S2n to the perturbations in kinetic parameters. The parameters perturbed are in 
the X axis and sensitivity of MK** (MK**-n) to such perturbations is shown in the Y axis. Only the most sensitive parameters and 'Sig' are shown. 
List of all the parameters of SI and S2 subjected to perturbation are shown in Additional file 2: Table SI with their respective values. List of all the 
parameters of Sin and S2n subjected to perturbation are listed in Table 3 of the main text. 
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perturbations in kinetic parameters in M3K layer. On 
the contrary, in S2 (Figure 9B) the incoming signal 
encounters the positive feedback before negative feed- 
back. Here the changes in the M3K layer are suppressed 
at the M2K layer by the negative feedback but as small 
changes in the MIC* can affect the strength of the posi- 
tive feedback at the M3K layer, the output MIC* exhib- 
ited maximum relative sensitivity to small changes in the 
MK layer itself (Figure 9B). S2n having identical archi- 
tecture of feedback loops as S2 also exhibited maximum 
sensitivity to changes in the MK layer and the layers 
below MK specifically to the shuttling rate of MK** be- 
tween the nucleus and cytoplasm (kllf and kllb in 
Figure 9D). 

Discussions 

Computationally it was predicted more than a decade 
earlier that MAPK cascade can exhibit oscillations em- 
bracing one negative feedback loop from MK** to sup- 
press M3K phosphorylation [14], much earlier than the 
experimental report on biochemical oscillations of the 
MAPK cascade [11,43]. Experiments have now shown 
that phosphorylation dynamics of MAPK exhibit oscilla- 
tory behavior from yeast to mammal [11,17,43]. Here we 
have studied the significance of differential designs of 
coupled positive and negative feedback loops in trigger- 
ing MAPK oscillations. We have also investigated how 
MAPK cascades embedded in designs such as PN-I and 
PN-II can shape their oscillation and the effect of 
nuclear-cytoplasmic shuttling of the cascade compo- 
nents triggered by each of the design. 

Oscillations in MAPK cascade due to PN-I and PN-II 
designs 

Although a single negative feedback is the minimal re- 
quirement for triggering MAPK oscillations, a growing 
number of studies indicates that oscillations in various 
cellular signaling systems [6,8] including the MAPK cas- 
cade [5], are triggered by coupled positive and negative 
feedback loops. These experimental reports led us to in- 
vestigate the roles of negative and positive feedback 
loops operative in a three-layer MAPK cascade (Table 1). 
Based on literature, we found that two possible designs 
of coupled positive and negative feedback loops can exist 
in a three layer MAPK cascade (SI and S2 in Figure 2), 
namely PN-I and PN-II. Our simulations show that both 
PN-I and PN-II can trigger oscillations in the cascade. In 
SI, the cascades output exhibited digital oscillations, 
whereas in S2 analogous oscillations were observed. 
These results show that the nature of the MK** output 
is determined by the type of the feedback loop func- 
tional in the M2K layer. From the context of information 
processing by a MAPK cascade, the ability to utilize two 
distinct designs of coupled positive and negative 



feedback loops would enable it to deliver unique oscilla- 
tory output while responding to input signal of similar 
strengths. We show that two MAPK cascades with iden- 
tical concentrations of their respective kinases and phos- 
phatases can trigger digital or analogous oscillations 
based on the design of coupled positive and negative 
feedback loop embedded in it. 

Information processing systems such as the signal trans- 
duction networks are usually activated by a spectrum of 
signals and strength of an incoming signal may not remain 
constant [41,44]. Thus in the living systems a signaling 
pathway needs to respond to signals of various strengths 
and subsequently deliver the desired output. We exam- 
ined whether the models SI and S2 can deliver oscillatory 
output when subjected to a wide range of signal strengths. 
It was found that both SI and S2 can exhibit their charac- 
teristic oscillations when subjected to a range of input sig- 
nal, although the system S2 was extremely robust to 
increase in signal strength above a threshold. The system 
SI exhibited equal amplitude oscillations whose oscillation 
frequencies were reciprocally dependent on the strength 
of the input signal (Figure 4A). However, S2 with feedback 
design PN-II exhibited equal amplitude and equal fre- 
quency oscillations for virtually any strength of input sig- 
nal, beyond threshold signal strength (Figure 4B). As the 
MAPK cascade is present in almost all the living systems, 
it is conceivable that the cascade is subjected to signal 
strengths varying in orders of magnitudes. We uncovered 
a remarkable ability of the cascade to trigger and maintain 
its oscillations with unchanged amplitudes and frequen- 
cies when subjected to varying signal strengths. A recent 
experimental report on epithelial cells stimulated with 
EGF also shows that the MAPK (ERK1) cascade conserves 
the frequency of oscillation of MK** when subjected to 
perturbations [27,28]. Our analysis reveals a plausible de- 
sign of coupled positive and negative feedback loops that 
the cascade can adopt to deliver such constant frequency 
oscillations. We additionally show that together with con- 
servation of amplitude, the cascade is also capable of pre- 
serving its oscillation frequencies in response to large 
fluctuations in incoming signals. 

Positive feedback emerging from an oscillating MAPK 
cascade triggers oscillations in its external target module. 

Literature of intra-modular crosstalk involving MAPK 
pathways is abundant. In T cell receptor triggered signal- 
ing pathways, MK** is the origin of 92% of the feedback 
loops (both positive and negative) [19], which implies that 
using the positive and negative feedback loops, the MAPK 
cascade determines fate of multiple pathways in the large 
scale network. Here we showed that oscillating MAPK 
cascade such as SI or S2 can use their respective positive 
feedback loops to trigger oscillations in any external signal 
transduction module. The extent of oscillation in the tar- 
get module would be determined by the ratio of rates of 
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phosphorylation and dephosphorylation in the target 
module. When the parametric conditions were satisfied in 
the target module (phosphorylation rate < dephosphoryla- 
tion rate), oscillations were triggered (Figure 5C). Oscilla- 
tions in the target module spanning from zero to its 
maximum phosphorylation amplitude were observed 
when phosphorylation rate was very much less than 
dephosphorylation rate. (Figure 5D). The ability to induce 
oscillations in the target modules depending on the ratio 
of kinetic parameters in the target module itself can be ex- 
tremely useful from the cellular context. This is because a 
plethora of target modules, each with unique ratios of 
phosphorylation-dephosphorylation will differentially de- 
liver their oscillatory outputs (Figure 5A-D). 

The result also exposes a multifaceted regulatory aspect 
of positive feedback loops which was not specifically 
addressed before. Positive feedbacks hallmark characteris- 
tics is signal amplification and promoting switch like be- 
havior to its target [25,33] . The feedbacks ability to trigger 
oscillations in its target (Figure 5A) reveals this novel 
regulatory aspect of the positive feedback. 

Fate of oscillations triggered by PN-I and PN-II upon 
nuclear cytoplasmic shuttling of MK layer and induction 
of its nuclear phosphatase 

Nuclear cytoplasmic shuttling of the MAPK cascade s MK 
layer components takes place and MIC* induces various 
transcription factors including its own phosphatases [10]. 
The models SI and S2 exhibited oscillations which are 
specific to cytoplasm but as MK layer of the cascade shut- 
tles between the nucleus and cytoplasm, fate of the oscilla- 
tions under such conditions is worth analyzing. We 
modified the oscillating systems where the modified sys- 
tems were built with both cytoplasmic and nuclear com- 
ponents (SI becomes Sin and S2 becomes S2n). The 
nuclear reactions comprised shuttling of MK, MIC and 
MIC* between cytoplasm and nucleus, P3-n induction fol- 
lowed by dephosphorylation of MIC*-n and MIC-n in the 
nucleus by P3-n. As the oscillations were triggered by the 
two different designs of feedback, PN-I and PN-II, we 
investigated how nuclear-cytoplasmic shuttling and tran- 
scriptional induction of P3-n affect the oscillations of Sin 
and S2n. Simulations show that oscillations triggered by 
the feedback design PN-I in Sin remains unaffected by 
the shuttling process and P3-n mediated dephopshoryla- 
tion in the nucleus (Figure 7A-D). However oscillations in 
S2n were abolished when nuclear phosphatase P3-n was 
transcribed in the nucleus (Figure 8A-D). Hence we show 
for the first time that fate of oscillations in a MAPK cas- 
cade is determined by the design of coupled positive and 
negative feedback loops that trigger such oscillations 
especially when compartmentalization of the cascade 
components take place. The study exposed probable cellu- 
lar strategies underlying generation and maintenance of 



robust MAPK oscillations for a longer duration, as long 
duration signal processing involves such nuclear cytoplas- 
mic shuttling and activation of various transcription 
factors. 

The feedback designs PN- and PN-II differentially 
determines the MAPK cascade's sensitivity to small 
perturbations in the model kinetic parameters 

Local sensitivity analysis was performed to understand 
the responses of the outputs MIC* (SI and S2) and 
MIC*-n (Sin and S2n) to small perturbations in their 
kinetic parameters (Figure 9A-D). Sensitivity analysis 
exposed the most sensitive parameters in the models 
embedded in the designs PN-I and PN-II. We found that 
sensitivity of MK** and MIC*-n exhibits differential sen- 
sitivity profiles in SI (Sin) and S2 (S2n), implying that 
the outputs sensitivity were determined by the design of 
the embedded feedback loops in the MAPK cascades. 
Sensitivity analysis results are useful for designing drugs. 
For example, for a system SI /Sin the most suitable 
strategy to suppress MK**/MK**-n will be to inhibit the 
strength of input stimuli (Sig) or enhance the flux of 
M3IC dephopshorylation. However if a drug needs to be 
designed for a MAPK cascade S2/S2n, MK**/MK**-n 
will be altered most effectively by altering the depho- 
sphorylation flux of the MK layer (for S2) or by altering 
the MK layer shuffling rates (for S2n). 

Proposed experimental verification of the model 
propositions 

The prediction made based on the simulation of the 
models SI, S2, Sin and S2n could be tested experimen- 
tally using different approaches. In the first approach 
mammalian cells such as COS-1 cells can be chosen to 
verify model type such SI. Experiments with COS-1 
show that MK** such as ERIC* gives positive feedback 
to M2K (MEK) phosphorylation step by inhibiting its 
competitive inhibitor RKIP [45]. At the same time ERIC* 
gives negative feedback to M3K (Raf) phosphorylation 
by inhibiting the upstream signal that triggers Raf phos- 
phorylation [5]. The design resembles the system design 
PN-I (SI) which also exhibited oscillations, as observed 
experimentally [5]. Hence considering COS-1 cells as ex- 
perimental system one could subject them with various 
perturbation conditions as described in the models. For 
example it is predicted from the simulations that SI can 
deliver oscillations with conserved amplitudes whose fre- 
quencies will vary according to the strength of incoming 
signal. Western blot analysis could subsequently be per- 
formed where kinetics of ERK phosphorylation for vari- 
ous strengths of input stimuli can be compared, which 
would then verify the model predictions. Further the 
model predicts that Sin (with design PN-I) should retain 
its oscillations upon nuclear-cytoplasmic shuttling and 
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induction of phosphatase such as MKP-1 should not 
affect the ERK oscillations. This can be tested by sub- 
jecting the COS-1 cells to prolonged stimuli and subse- 
quently capturing the phosphorylation kinetics of ERIC*, 
which should exhibit oscillations, as predicted by the 
simulations. Presence of oscillations during the nuclear 
cytoplasmic compartmentalization of the ERK cascade 
can be experimentally tested in the same lines as 
explained elsewhere [11]. 

The system design S2 where positive and negative 
feedbacks are coupled as design PN-II are not reported 
in one single study as yet. But a recent study shows that 
three layer MAPK cascade can be synthetically built 
[44]. Such synthetic systems will be ideal for testing hy- 
pothesis. One could design the system S2 as a synthetic 
system. Mass spectrometry data suggest that ERK** pro- 
vides positive feedback to Raf (M3K) by phosphorylating 
it in certain residues which enhances specificity of Raf 
phosphorylation by many fold [46]. Coupled to that a 
negative feedback from ERK to Raf can be considered in 
which ERK hyperphosphorylates and desensitizes Raf 
[47]. The overall design would resemble the system de- 
sign PN-II. Here the positive feedback is in the form of 
enhanced Raf phosphorylation in response to the incom- 
ing signal which is followed by the negative feedback in 
the form of desensitization of phosphorylated Raf 
(M3K*) that will consequently inhibit MEK (M2K) phos- 
phorylation. Such synthetic cascades with positive and 
negative feedback resembling design PN-II could be sub- 
jected to signals of variable strengths and the oscillatory 
amplitudes of the cascade output can be captured in the 
form of western blots. The simulations proposed that 
the system S2 subjected to a very wide range of input 
signal should exhibit oscillations with conserved ampli- 
tude and frequencies which could be verified building 
the synthetic MAPK cascade. 

Biological significance of MAPK oscillations and proposed 
implications of this study 

Exact biological message encoded in the oscillatory 
waves of the MAPK cascade is not yet understood well, 
though it is argued that the oscillatory MIC* fulfils some 
requirement for triggering transcription of certain cyclic 
genes [48]. The current archetype states that, signaling 
system in general encodes messages either in amplitude 
or in frequency (or may be in both) of the oscillatory sig- 
nals, for triggering transcription of a plethora of genes 
[11,12,24,48]. Here, through our study we demonstrated 
various ways in which unique oscillatory message could 
be transmitted by the MAPK cascade embedded in 
coupled positive and negative feedback loops to its nu- 
clear targets. 

The feedback design PN-I can trigger oscillations of 
equal amplitudes but of variable frequencies. This type 



of cascade could be utilized by the cell for activating a 
subset of cyclic target genes, all of which require identi- 
cal amplitude of MK** as their activation threshold but 
the interval of expression of each target gene is deter- 
mined by the frequency of oscillations. The feedback de- 
sign PN-II can be utilized to deliver oscillations with 
near identical frequency and amplitude in response to 
widely varying signal strengths. This type of feedback de- 
sign would be suitable for a MAPK cascade involved in 
robustly inducing specific sets of genes whose expres- 
sions are critically dependent on the amplitude and/or 
frequencies of the MK**. 

We demonstrated how oscillations could be main- 
tained during a long duration signaling when signal pro- 
cessing involves nuclear-cytoplasmic shuttling of the MK 
layer of the cascade, followed by transcriptionally indu- 
cing the phosphatases that interact with the cascade it- 
self. We showed that it is not always possible to 
maintain oscillations in the face of obvious biological 
perturbations, such as interaction with the transcription- 
ally induced phosphatases and thus the cascade has to 
adopt certain feedback designs (such as PN-I) to endure 
such perturbations to exhibit prolonged oscillations. 

Conclusion 

The MAPK cascade can utilize architecturally distinct 
organizations of coupled positive and negative feedback 
loops to trigger its oscillations. We uncovered that the sig- 
naling pathways such as the MAPK pathway can uniquely 
process wide range of signals by utilizing its feedback 
loops. It is intriguing how adoption of specific design (PN-II) 
of coupled feedback loops can trigger oscillations with ex- 
tremely robust frequency and amplitude, specifically when 
such robustness in the oscillations are desired in an envir- 
onment where the external signal strength fluctuates in 
several orders of magnitudes. Subsequently we show the 
trade off associated with such feedback designs (PN-II) 
during the nuclear cytoplasmic compartmentalization of 
the cascade, where oscillations triggered by PN-II couldn't 
sustain such compartmentalization effect. However oscilla- 
tions triggered by PN-I were robustly maintained during 
the compartmentalization of the MAPK cascade compo- 
nents. Thus it can be argued based on our analysis that 
MAPK cascade embedded in PN-II can be used by specific 
cell types to exhibit short duration oscillations in response 
to extremely noisy signal, where frequency and amplitude 
needs to be robustly maintained. The oscillations triggered 
by PN-II will be of short duration as longer duration in sig- 
naling implies nuclear compartmentalization of the MAPK 
cascade, which leads to attenuation of PN-II triggered 
oscillations. On the contrary the design PN-I can trigger 
long duration oscillations (involving nuclear cytoplasmic 
compartmentalization), when the cascade embedded in 
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such design is exposed to a relatively less noisy input 
signal 

We additionally found a completely unexpected regu- 
latory behavior of the positive feedback component of a 
coupled positive and negative feedback loop used for 
triggering MAPK oscillations. We show that positive 
feedback emerging from an oscillating MAPK cascade 
can generate a spectrum of unique oscillatory informa- 
tion to various external target modules. The amplitude 
of oscillations thus triggered would depend on the ratio 
of phosphorylation and dephosphorylation in each of the 
target modules, which means, each target can attain dif- 
ferential oscillatory fates by adjusting such ratios. 

Availability of supporting data 

The supporting data are provided as additional files with 
the manuscript. Additional files include Additional file 3: 
Figure SI, Additional file 3: Figure S2, two Additional 
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induction of P3-n was stopped. P3 was present in the system (P3 = 500 nM) 
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